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A thin layer of liquid in a horizontal cell is subjected to a periodic vertical force with 
two control parameters: acceleration and frequency. The influence of the rheological 
behavior of the fluid was considered over the empirically obtained results, which 
were subjected to a multfractal spectrum process. As control parameters varied, 
so did the corresponding spectra, according to theoretical models; each change of 
state was visually detected, which permitted interpreting the corresponding spec- 
tral changes. 
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1. Introduction 

The classical Faraday patterns are a selection of surface waves following a spe- 
cific criterion studied by Benjamin and Ursell (1954). The amplitude and the wave 
number of the surface follow a certain harmonic oscillator equation. In this paper, 
the liquid is shear thinning, adding further complexity to the classical case. A thin 
layer of liquid in a horizontal cell is subject to a periodic vertical force, its frequency 
and acceleration being the control parameters. Between the acceleration/frequency 
threshold (acceleration measured in g units — g being the acceleration due to grav- 
ity) after which the classical Faraday patterns appear, and the new threshold cor- 
responding to the ejection of droplets (spray), there is an intermediate zone with 
a structure that is the object of our study. This intermediate state is characterised 
by the formation of 'cusps' on the surface of the liquid under study; the behaviour 
of the geometrical patterns formed by these 'cusps' depends on the control param- 
eters of the system. The liquid used in these experiments is a polymeric solution of 
polyacrylamide (PAAm) in a 2g/l concentration. PAAm is a shear-thinning poly- 
mer, i.e. there is a region where its viscosity diminishes with the shear rate. In 
this paper, the behaviour of the 'cusps' is studied by imposing a voltage difference 
across a pair of electrodes. One electrode is submerged in the liquid; the other is 
suspended over the liquid surface at a controllable distance. When a cusp rises and 
touches the electrode above, it closes a circuit and a current is generated yielding a 
signal, electronically detected and processed. In figure 1 we present the experimen- 
tal setup and figure 2 shows the images corresponding to three different frequencies 
and the response of the detector, where we can observe the signal as a function of 
time. With constant frequency v, we vary the acceleration n and obtain a signal as 
a function of time. We study C K , the Cantordust intersection of the signal with a 
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segment. 

In the literature there are two ways to obtain a multifractal spectrum (al- 
pha, f( alpha)) of a Cantordust: one is through the very definitions of the alpha- 
concentration and the corresponding spectral value f (alpha), the other is through 
the use of the lagrangian coordinates. It is not proved that the two methods yield 
the same spectrum. The spectrum given by the lagrangian coordinates is an ideal 
one, to which the spectrum given by the definitions is supposed to approximate. 
Yet, the latter can have deviations from the ideal spectrum. In this paper, we work 
with the definitions (despite the difficulties involved, e.g. far less points in the spec- 
tral curve), since such deviations from the ideal shape give useful information on 
the nature of the experiment. 



2. Experimental setup 

The experimental setup consists of a square acrylic container (side 7 cm), on top 
of a vibrator. The system is excited in the vertical direction with an acceleration 
k = ao sin vt, where v is the driving frequency and ao the other control parameter. 
A pair of electrodes, one submerged at the bottom of the cell, the other 3 mm from 
the surface, permits detecting the contact of the liquid with the second electrode, 
thus closing an electric circuit and sending a time signal. 

The frequency applied varies between 20 and 120 Hz, the thickness of the liquid 
layer being 4 mm. 

As noted in the Introduction, the polymeric solution is a shear thinning liquid. 
We can observe (Fig. 3) that the viscosity diminishes with the shear rate. In Fig. 2 
we observe the formation of cusps, and intermediate stage between pattern forma- 
tion and droplet ejection. 

With constant frequency v, we vary the acceleration k and obtain a time signal. 
C K is the Cantordust intersection of the signal with a segment. 
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Figure 2. Cusp images and obtained signal 
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Figure 3. Viscosity (centipoise) as a function of shear rate (1/s) for the experimental 

liquid PAAm, 2g/l 



3. Multifractality in an empirical context 

We will work in empirical fractal sets F contained in the unit segment [0, 1], en- 
dowed with a probability measure \x. Let us recall that the a-concentration of a 
box B C [0, 1] of length I is the log/log version of the density — mass or measure 
per unit length — that is a(B) — log ^J ; nB . The a-concentration of a point in F 
is defined (just like the density of a point) as the limit of a's of shrinking boxes 
containing said point. The spectrum /(a) is the box dimension of the elements in 
F that share the a-concentration. 

When dealing with an experimental sample, an empirical Cantordust fractal 
set embedded in a unit segment, which we want to process in order to obtain a 
multifractal spectrum (a, /(a)), some practical difficulties appear. There are sev- 
eral rather severe discrepancies with the theory proper. On the one hand, a fractal 
Cantordust-like set has a non-numerable infinity of elements, whereas in any con- 
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Crete empirical context, we have a finite number of measurements in the sample. 
Hence, the theoretical concept "e — ► 0" is rendered meaningless: say that S is the 
size of the sample taken at regular intervals, and recall that the Cantordust is em- 
bedded in a unit segment. Then ^ = As = e s , where "s" stands for "sample", is 
the smallest possible value of such e. Next, we cover the unit segment with boxes of 
a certain length Al in order to define the so-called natural probability measure /i. 
Theoretically, such natural measure is defined as the proportion of sample elements 
in a box of length Al = ei "as e; — > 0" . . . again impossible in an empirical context. 
In such a context we must have e s considerably smaller than e;, i.e. e s <C ej. There- 
fore we need to have S>B, the total number of boxes of equal size. 

There is a yet more delicate e to take into account: having e s and ei, e s <C e/, and 
the natural measure u, we establish the values a for each box B,a = l ° s ^ BnF ^ . 
We obtain thus a number of a- values, from a certain a m i n to an a max . Theoreti- 
cally, we have to consider f(a) = lim ei ^ *± , N a the number of elements in the 

infinitesimal segment [a, a + e Q ), e a arbitrarily small and tending to zero. It may 
well be that the (finite number of) a- values empirically obtained between a m in and 
ctmax, are all different. . . in which case, if Aa = e a is rather small, each interval 
[a, a + Aa) may have one element (if any). . . in which case /(a) = 0. 

In practice, it should be £ s < t; < e„, though in theory all these e tend to zero 
independently of one another. This implies that S » B which in turn should be 
3> A, the number of points (a^, /(««)) obtained in order to form a smooth curve 
(a,f(a)). We have found, elsewhere (Piacquadio and Cesaratto 2001) and in this 
work, that in this context "3>" means "larger than the square of, e.g. S should 
be larger than the A th power of A. For instance, when calculating the multifractal 
spectrum of the Farey tree (Piacquadio and Cesaratto 2001) we needed an S of the 
order of 10 6 , whose A th root should yield 31 points (c^, /(«»))• Yet, we could obtain 
only some 15 such points constituting a smooth curve (a, f(a)) that agreed with 
the theoretical spectrum. 

The empirical context of this work admits a maximum value of 10 4 elements in 
S, for beyond that magnitude the experiment changes its nature, as seen in the last 
section. 

Given that S = 10 4 , then VW = 10 2 has to be the maximum number B of 
boxes — the absolute maximum, a slightly smal ler n umber would be more adequate. 
Hence, a number slightly smaller than \/B = VlO 2 = 10 points (a*, /(o^)) will yield 
a smooth curve (a, /(a)) in each of the cases below. 



4. The evolution of the Cantordust set C K 

In order to analyse the evolution of C K as k increases, we adapt the model of 
the Chaos diagram for the one-parameter Henon-like system (Grebogi et al. 1987; 
Alligood and Sauer 1993; Robert et al. 1999). In this diagram the orbits of a one- 
parameter — denoted here by k — family of Henon-like maps are represented, by 
means of the global bifurcations they undergo, as n slowly increases. Chaos is 
reached for a certain specific value k$, which means that the system is charac- 
terised by the presence of a one-piece connected fractal attractor F K . As variable 
k slowly increases beyond Ko, the system in chaos undergoes certain key changes. 
Grebogi et al. (1987), Alligood and Sauer (1993) and Robert et al. (1999) studied 
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Figure 4. The spectrum f K {a) for k < k c = 4.93(7 an d frequency v = 20Hz 

the dynamics of the change — the evolution of the so-called accessible rotation num- 
bers for the one-piece fractal attractor F K . When the increasing variable k reaches a 
certain critical value k c , the chaotic system breaks down, and undergoes a so-called 
boundary crisis, in which the attractor F K suddenly disappears. Some work has 
been done (Hansen and Piacquadio Losada 2005) in order to relate the evolution 
of the dynamics of this chaos-on-road-to-crisis model to the geometrical changes 
in the one-piece F K , as k approaches k c , in which F K disappears and the system 
breaks down. Following this lead, we analyse the evolution of C K by means of that 
of the geometry of f K = f K (a), the multifractal spectrum of the natural measure 
[i K associated with C K . The guiding principle may be stated as follows: we study 
fractal C K by means of its natural measure /^ K , [i K by means of its multifractal 
spectrum f K , and the evolution of the system by the geometrical evolution of the 
curve (a, /«(«)). 

Briefly, details follow, the chaos-on-road-to-crisis type of situtation, as k nears 
a critical value k c , is characterised by a one-piece connected multifractal spectrum 
f K , which fulfills all the right thermodynamic laws. When k reaches the value k c the 
system reaches a crisis point and breaks down, as the Lagrangian coordinates fail to 
make sense. For k > n c the measure ji K seems to disappear, for the corresponding 
multifractal algorithm produces a function made up of unconnected pieces of what 
looks like fragments of different spectra — so the thermodynamic laws cease to be 
valid. 

5. Chaos on road to Crisis 

For k near k c — 4.93g, e.g. for n = 4.08.g, the natural measure \i K of C K produces 
a one-piece connected multifractal spectrum, as seen in Fig. 4. C K was produced 
by a sample of size S = 10 4 ; B = VS — 100, and A slightly smaller than a/B- We 
observe that the spectrum is convex, that there are points in the graph at both left 
and right of a = 1, that the maximal dimension is rather high — above 0.7 — in 
agreement with the naked-eye observation of the sample. The one law which does 
not appear to be fulfilled is f(a) being tangent to the line y = x, as shown in Fig. 
5. This is due to the small value of B. We increase the number B of boxes — setting 
at risk the smooth character of the graph of the spectrum — to 200. We observe 
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Figure 5. The spectrum f K {a), k — 4.08g, 100 boxes. 
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Figure 6. The spectrum f K (a), k — 4.08g, 200 boxes. 
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0.7457 
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1.1278 



Table 1. The spectrum f K (a), k = 4.08<7, 100 and 200 boxes 

that the convex character of the curve is preserved, as seen in Fig. 6, together 
with the features described above for B = 100. In table 1 we compare a m i n for 
B = 100 and B = 200, we obtain a mm (200) = 0.87 • • • < a mm (100) = 0.90 • • • . The 
same happens for the values of au for which the spectrum reaches its maximum: 
a M (200) = 0.94- •• < a M (100) = 0.96 ••• . Both values indicate that, when B 
increases, the spectrum shifts towards the left, approaching the bisectrix line y = x. 
Also, .f maa; (200) = 0.74 • • • > / max (100) = 0.72 • • • , which indicates that the shift is 
upwards as well as leftwards, getting sensibly closer to y = x. Increasing B further 
already distorts the smooth character of the spectrum. 
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Figure 7. The spectrum / re (a) for k = k c = 4.93g and frequency v = 20Hz 
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Figure 8. The spectrum / K (a), k = n c — 4.93g, 100 boxes. 



6. Crisis: critical value k c of k 

The frequency v = 20 Hz is kept constant. When k reaches the critical value 
k = k c — 4.93g, keeping S = 10 4 and B = vS = 10 2 , we obtain for f(a) = 
fn c ( a ) a graph with f'(a) constant along a large segment, as shown in Fig. 7. 
The segmented part of the graph implies a breakdown of the thermodynamical 
formalism, for the Lagrangian coordinate q = f'(a) collapses: if q is constant, 
it cannot be the parametric coordinate of both a and f(a) shown in the figure. 
The segmented characteristic corresponding to this value k c of k is stable under 
variations of B: we observe (see Fig. 8) that, as before, the graph is not tangent to 
the bisectrix y = x. We increase B to 150 (sec Fig. 9), some irregularities appear 
in the curve; still the segmented part remains quite apparent. Although the values 
of OL m in and a max for the two graphs of f Kc (ct), for B = 100 and B = 150, are 
much the same, as seen in Fig. 10, Table 2 shows that for B = 150 the curve moves 
upwards, nearing line y = x. Moreover, ajvf(B = 100) > cva/(B = 150), which 
implies that the upper cusp of the curve moves towards the left (when B grows), 
again nearing y = x. As before, the "distance" between the spectrum and line y = x 
is only apparent, due to the low value of B. For B = 200 we still have the segment, 
and the graph approaches even more the bisectrix,. . . but the curve begins to look 
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Figure 9. The spectrum f K (a), K — n c — 4.93gr, 150 boxes. 
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Figure 10. The spectra f K (a), for k = k c = 4.93g, 100 and 150 boxes. 
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B = 150 



OtKiifmax) 
/W = fmax 



0.9788 
0.7235 



0.9762 
0.7923 



Table 2. Maximal values for f Kc (a), B = 100 and 150. 



shaky, losing its smoothness. 

The breakdown of the system, as source of information, goes beyond the collapse 
of the Lagrangian coordinate — i.e. /'(a) — constant. The range of internal energies 
a allowed is Aa = ct max — a m i n = 1.133 — 0.943 = 0.19, whereas from Table 1, Aq 
is larger than 0.22: the range of allowed internal energies is minimal for n = k c , 
restricting further the information obtained from the spectrum. Also, ajvr, for which 
f(a) reaches its maximal value, is much nearer unity for k = n c than for k — 4.08<7, 
as seen in Tables 1 and 2. The value a = 1 yields no information, since it means a 
uniform distribution of probability vis-a-vis length, hence, a minimal variation in 
such distribution. 
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7. Segmentation and the Henon map 

A multifractal spectrum with a possible segmented part was studied by Gunaratne 
and Procaccia (1987): the Henon map. In this system, anomalous periodic orbits 
with atypical behaviour contribute anomalous and isolated values of a, and the end 
of the smooth part of the f(a) spectrum can be connected to some isolated (a, f(a)) 
point by a segment; /'(a) remaining continuous. The anomalous cycles contribute 
to isolated measures which, in the global partition sum from the thermodynamical 
formalism, yield such isolated (a, f(a)) that appear to suggest a segmented portion 
of the global spectrum. In the Faraday experiment we find that f Kc (a) is segmented 
because other measures get integrated in the global partition: for k > k c we find 
more than one measure clearly identified, contributing to the global multifractal 
algorithm. 

In the Henon case, the segmented f(a) is conjectured to represent a possible 
phase transition — as previously suggested by Grassberger, Badii, and Politi (pers. 
comm. to Gunaratne and Procaccia). This segmented situation is deemed to be 
quite marginal, being extremely sensitive to minute variations of the Henon map's 
parameters. We interpret the segmented curve in Fig. 7 much in the same terms, 
for f(a) = f Kc (a) is a marginal state extremely sensitive to minute variations of 
k, since for k < k c there is no segment, and for k > k c there is no /i^, but two 
distinct measures (bi-multifractality, as seen below) which can be interpreted as a 
phase transition, extending the conjecture of Grassberger, Badii, and Politi. 

8. Post-crisis case k > k c : Bi-multifractality 

Radons and Stoop (1996) studied the phenomenon of bi-multifractality obtained 
when two measures compete in the same support. They superposed two indepen- 
dent strictly self-similar multifractal measures supported in two Cantordust sets in 
the same segment. The two measures yield independent spectra fi(a) and /2(a), 
theoretically known, since they are self-similar. The references quoted in their work 
relate previous attempts at such superpositions. In Fig. 3b) in their paper, we ob- 
serve the three spectra: the theoretically known fi(a) and /2(a) in faded dotted 
lines and, in solid lines, the global generalised thermodynamical formalism that 
yields the empirical f(a), obtained by integrating the data from both measures: 
covering the whole segment that supports both measures with boxes and carry- 
ing out the multifractal algorithm. The resulting global integrated f(a) traces a 
bit of (a,fi(a)) and a bit of (a, /2(a)), with a substantial gap between the two 
broken, unconnected fragments: bi-multifractality. The graphs of fi(a) and /2(a) 
intersect at two points Pi and P r , "Z" and "r" for "left" and "right" . One of the 
stopping points of the fragmented global f(a) is P r , where /{ and f 2 are finite and 
have different signs. This is interpreted by Radons and Stoop as a change of phase, 
second-order transition. Fig. 11 shows global integrated spectra taken from the 
Faraday experiment, for k > n c in each case. Fig. 11 a) shows an embryonic /2(a). 
Figs, lib) and c) show unconnected fragments of two multifractal spectra, sepa- 
rated by a gap: bi-multifractality. Neither of the two fragments of spectrum shows 
the symmetries corresponding to self-similar multifractal characteristics. But the 
geometrical features common to the pairs of fragments in Figs. 11 a), b) and c), 
suggest a sketching of a bi-multifractal configuration — see Fig. 12 — in accordance 
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Figure 11. Fragmented spectra (bi-multifractality) for high values of re 




Figure 12. Bi-multifractal configurations — Figs. 11(b) and 11(c) — in accordance with 

the Radons and Stoop model 



with the second-order transition of the Radons and Stoop model. The shape of 
the fragmented spectra f K ,K>K c in Figs, lib) and c), corresponding to different 
cases, suggests a common road post-crisis-through-bi-multifractality for this Fara- 
day experiment. 

In the Henon case, the fractal attractor F K disappears for k > k c ; analogously, 
the post-crisis type of situation for k > k c in the Faraday experiment shows that the 
multifractal measure /i K disappears — for the set of thermodynamical laws ceases to 
be valid for the broken f K (a). This disappearance of fi K is explained by the presence 
of two (or more?) conflicting measures. 

There is a continuity among the chaos-road-to-crisis, crisis, and post-crisis-bi- 
multifractality cases: the first one, k < k c , shows a "normal" fi K through a "normal" 
/(a); the second shows anomalous measures, yielding a not normal 

spectrum through the appearance of a segment, which is interpreted as a change 
of phase. The third case shows, explicitly and clearly, those conflicting measures: 
the spectrum is not "anomalous" any more, but broken down, the phase transition 
now quite apparent. 
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9. Qualitative observations 

Past the crisis point, we obtain, in the three cases of Fig. 11, a bi-multifractal 
situation. The three cases seem to share some qualitative — if not quantitative — 
properties. The case of Fig. 11 a) is the weakest, since the second multifractal 
reduces itself to a single point. Yet, it is a point that cannot be ignored: on the one 
hand, it is not on the real a— axis, which means that it is not on a freak isolated 
box with just one lonely clement in the Cantorset, the other 9,999 being together 
elsewhere. On the other hand, it is far enough away from the typical and usual 
(a, /(a)) at its left, as to make it impossible to consider this point as part of that 
convex curve. The one-point spectrum is exactly at the height of one of the points in 
the typical graph at the left, a feature also present in the other two graphs — Figs, 
lib) and c). With higher values of the parameters, the two broken spectra (Fig. 
lib)) appear more clearly, yielding bi-multifractality. Fig. 11 c) presents the same 
type of bi-multifractality (from a qualitative point of view), for the same liquid 
at a higher viscosity. A word on the significance of this fact: as an example, let 
us consider the Farey tree (Piacquadio and Rosen 1995, see also Piacquadio and 
Cesaratto 2001). The tree is a partition of unit segment [0, 1] — [f, y], obtained 
by interpolating rational number between | and ^7. Successive interpolations 
give rise to all rational numbers, starting from \ = y±i interpolated between 2 and 
j. A certain optoelectronic experiment (Essevaz-Roulct et al. 2000) shows a camera 
that films a screen to which it is connected, filming its own image, feeding back 
the image to the screen. When the camera turns around its optical axis, a pattern 
of p light spots appear on the screen, after q turns. Integers p and q follow the 
hierarchy of the Farey tree. . . But, when measuring p's and q's, the precision of such 
measurements has to increase with q, so the Farey tree has a natural bound for the 
q's in each refinement of the partition: all rational numbers | will appear, but they 
do so slower than in the theoretical interpolation, for they do it, perforce, according 
to the size of q. All other things remain equal: the interpolation laws in the unit 
segment, the equidistribution of the probability measure. . . Yet, when processing 
the theoretical partition (Farey-Brocot tree) as a multifractal spectrum fp-B{oi), 
and the practical one /f(«), perforce a slower production of rational numbers — due 
to inescapable experimental reasons — we find that /f-b(oO and /f(cc) cannot be 
more different. The theoretical /f-b(cO is strictly increasing and the range of the 
a's is [a min ,oo): an infinite segment; whereas /f(«) is strictly decreasing, the a 
range is [1,2]. Hence, the same interpolation rules, the same measure distribution 
rules, applied to the same fractal support,. . . with an inescapable slow-down in 
the velocity of interpolation, yield totally unrelated multifractal spectra. Our point 
being that the algorithm yielding a multifractal spectrum is extremely sensitive to 
(apparently innocuous) variations in its measure. Hence, we can conjecture that, 
when two spectra have similar qualitative characteristics, as is the case illustrated 
in Figs. 11 b) and c), the underlying phenomena studied with multifractal tools 
have to be strongly connected. 



10. Discussion and Conclusions 



For different values of k, the processed data of the corresponding signals — accord- 
ing to the frame described above — agree with direct observations of the experiment: 
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Figure 13. Different viscosity/shear-rate curves imply that the liquid has become 
non-homogeneous for S > 10 4 . 

for k /* k c (chaos-road-to-crisis) we observe the cusps maintaining the same locus 
on the surface and the same height. For k > k c (bi-multifractality) the cusps are 
sparse and do not preserve their location on the surface. They group in clusters that 
'wander around' on the liquid surface. When we arrive at 30 Hz and k » k c the 
cusps are sporadic and appear to be random. Beyond bi-multifractality we observe 
the ejection of droplets (spray). For S > 10 4 the liquid starts consistently losing 
cusps and the system tends to change its structure, as shown in Fig. 13. 

The laboratory experiment perforce limits the size of the data set available for 
analysis. Despite this limitation, multifractal data processing permits the identifi- 
cation of distinct states (chaos, crisis, . . . ). These (laboratory) results are in accor- 
dance with pre-existing theoretical models of such states. Especially noteworthy is 
that the different states can be visually observed in the course of the experiment. 



We are indebted to Carlos Lizarralde for his technical assistance in measurements 
and data processing, for helpful suggestions, and for encouragements too numerous 
to mention. 
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